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Abstract 

We consider a size-structured aggregation and growth model of phytoplankton community proposed by Ackleh and Fitz¬ 
patrick Q. The model accounts for basic biological phenomena in phytoplankton community such as growth, gravitational 
sedimentation, predation by zooplankton, fecundity, and aggregation. Our primary goal in this paper is to investigate 
the long-term behavior of the proposed aggregation and growth model. Particularly, using the well-known principle of 
linearized stability and semigroup compactness arguments, we provide sufficient conditions for local exponential asymp¬ 
totic stability of zero solution as well as sufficient conditions for instability. We express these conditions in the form of 
an easy to compute characteristic function, which depends on the functional relationship between growth, sedimentation 
and fecundity. Our results can be used to predict long-term phytoplankton dynamics 

Keywords: Nonlinear evolution equations, principle of linearized stability, spectral analysis, 
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1. Introduction 


Planktonic lifeforms provide a crucial source of food (organic carbon) for many aquatic species including blue whales 
E3- In fact, oceanic plankton collectively provide approximately 40% of worlds organic carbon |8). The main component 
of plankton community are unicellular algae called phytoplankton flil ]. Similar to terrestrial plants, phytoplankton make 
their living by photosynthesis, and consequently inhabit surfaces of lakes and oceans. Phytoplankton cannot swim against 
a current and thus form aggregated communities at the surfaces of lakes and oceans to promote survival and proliferation. 
Besides their macroscopic predators such as whales and shrimp, aggregated phytoplankton community are also removed 
due to gravitational sedimentation and zooplankton grazing on phytoplankton. 

To study the dynamic nature of phytoplankton communities various mathematical population models have been 
developed. The mathematical model that we consider in this article is the model first considered by Ackleh and Fitzpatrick 
in 2]. Although, the existence and uniqueness of a global positive solution in several different spaces have been proved in 
MS, the long-term behavior of this model has not been investigated. This is mainly due to the nonlinear nature of 
the Smoluchowski coagulation equations used for modeling aggregation. Hence, our main goal in this paper is to derive 
sufficient conditions for the local exponential asymptotic stability of the zero solution for the aggregation-growth model 
(see Section 0. In Section [3l we also derive sufficient conditions for instability of the zero solution. These conditions can 
then be used to predict long-term phytoplankton dynamics. 

The size-structured population model proposed in [2] models growth of aggregates due to cell division and aggregation 
and removal of aggregates due to sedimentation and microscopic predation by zooplankton. In a phytoplankton community, 
the density of aggregates of size x at time t is denoted by p(t, x). An aggregate is assumed to have minimum xq and 
maximum x\ possible sizes. In vivo, there are no aggregates of volume 0 and the aggregates cannot grow indefinitely, so 
the only biologically plausible case is0<a;o<a;i<oo. Hence, in this paper we consider the case 0<*o < x\ <oo, and 
postpone the analysis of the case with xq = 0 and/or xi = oo for our future papers. 

We will consider the following nonlinear partial integro-differential equation model for the evolution of a phytoplankton 
population, 

d t p(t, x) = - d x (gp ) + F\p], g(x 0 )p(t, x 0 ) = K[p\(t), p{ 0 , x) = p 0 (x) £ L^xq, x\] ( 1 ) 


where 


and 


^ px—XQ r-X\ 

r\p\{t, x) = - P(x-y, y)p(t , x-y)p{t , y)dy-p(t , x ) / /3(x, y)p(t, y) dy - w(x)p(t, x) 

^ J Xq J Xq 

£[p](*) = [ q(x)p(t, x)dx. (2) 

J Xn 
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The function g(x ) represents the average growth rate of the aggregate of size x due to mitosis. Specifically, when a single 
cell in the aggregate of size x divides into two identical parent and daughter cells, the daughter cell enters the aggregate of 
size x contributing in a increase in total size. The coefficient w{x) represents a size-dependent removal rate. Biologically, 
aggregates can be either removed by gravitational sedimentation or zooplankton grazing on phytoplankton mi- /?(x, y) 
is the aggregation kernel, which describes the rate with which the aggregates of size x and y agglomerate to form an 
aggregate of size x + y. The fecundity rate q{x) in @ represents the number of new cells that fall off an aggregate of size 
x and enter single cell population. 

As our solution space we use H = T 1 [a"o, Xi] with the usual norm ||-|| L i (hereafter, just ||-||). Consequently, the 
equation 0 can be written as a semilinear abstract Cauchy problem (ACP) 


Pt = C\p] + A f\p], p( 0, x) = bo{x) G H. 

The operator C : T>(C) C H —¥ H is defined as 

C[p\{x) = - (g(x)p(x))' - w(x)p(x) 

with its corresponding domain 

V(C) = {(/>£ H\ (#)' G H, (#)(x 0 ) = JCm . 

The nonlinear operator M : H —H is defined as 

^ nX — Xo rX i 

Af[p] = - / P{x-y, y)p(x - y)p(y) dy - p(x) / /3(x, y)p{y) dy. 

^ J Xo J Xq 

We make the following assumptions 

g G C' 1 [*o, Xi]; and g(x) >0 for xq < x < x\, 

/3 G L°°([x 0 , xi] x [x 0 , xi]); /3(x, y) = f3(y, x) and f3(x, y) = 0 if x + y > x x 

w G C[x o, X\] and ui > 0 a.e. on [xo, X \], 
q G L°°[x o, x{\ and q > 0 a.e. on [xq, Xi] . 


(3) 

(4) 

(5) 

( 6 ) 


(•*4) 

K/?) 

(^ 4 ) 

(•<) 


Note that the restriction on g states that any aggregate of size x G [xo, Xi] has strictly positive growth rate. Any aggregate 
growing out of the bounds are not considered in the model. The assumption on w(x) enforces continuous dependence of 
the removal on the size of an aggregate and ensures that any aggregate size has a non-negative removal rate. Having the 
required ingredients in hand, we present the main result of this paper below, and demonstrate our proof in the subsequent 
sections. 


Theorem 1. Under the assumptions (■srfg), (&&w) and (srfq), the zero solution of the nonlinear evolution equation 

defined in 0 is locally asymptotically stable if 



Moreover, the zero solution is unstable if 




dx < 1. 


dx > 1. 


2. Linearized stability for the zero solution 


Our main goal in this section is to give sufficient conditions for the zero solution, p[t , x) = 0, to be locally asymptotically 
stable. In particular, we use the linearized stability results introduced in 0- For the convenience of readers we now 
present several results, which can be found in most of the semigroup theory books (see Q for reference). 

The growth bound, wo(A), of a strongly continuous semigroup (T(t)) t>0 with an infinitesimal generator A is defined 


as 


ujq(A) := inf < w G 


3./14 > 1 such that 
||T(f)|| < M u e ult for all t > 0 


In the development below we denote DA(f) to be the Frechet derivative of an operator A evaluated at /, which is defined 
as 


DA{u)h = A[u + h] — A[u ] + o(h), Vit G T>(A ), 


where o is little-o operator satisfying ||o(/i)|| < b(i ■) \\h\\ with increasing continuous function b : [0, oo) —> [0, oo), 6(0) = 0. 

The stability results of this section are based on the following proposition from 0 . p. 198]. We refer readers to [ 3 ] for 
a generalized version of this proposition, which applies to a broader range of nonlinear evolution equations. 
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Proposition 2. Let (T(t )) t>0 be a Co semigroup in the Banach space X with an infinitesimal generator C. Let M : X —> 
X be continuously Frechet differentiable on X. Let f € P(£) be a stationary solution of ([5ji. i.e., (C + M) [/] = 0. If 
the linearized operator C + DJ\f(f) (which is the infinitesimal generator of a Co semigroup by the well-known perturbation 
theorem) satisfies u>o (£ + LW(/)) < 0, then f is locally asymptotically stable in the following sense: 

There exists rj, C > 1, and a > 0 such that if ||/ — /|| < y, then a unique mild solution u(t) of ([5]), 

u(t) = T(t)f+ f T(t — s)J\f[u(s)] ds , 

Jo 

exists for all t > 0 and ||w(t) — /1| < Ce~ at ||/ — /|| for all t > 0. 

In jij authors have proved that the linear operator C generates a strongly continuous semigroup. Hence, we prove the 
second assumption of Proposition [2] below. 

Lemma 3. The nonlinear operator J\f defined in Q is continuously Frechet differentiable on H. 

Proof. The Frechet derivative of the nonlinear operator J\f is given explicitly as 


2 pX — Xq nX\ nX\ 

DN{4>)[h{x)) = - / P(x-y,y)[<j>{y)h{x-y) + h(y)cf>{x-y)\dy-h(x) j /3(t, y)<^(?/)ciy-</>(x) / /3(x, y)h(y)dy. 

^ J X 0 J Xq J Xq 

For the arbitrary functions u\, U2 G H we have 

\DJ\f(ui)h(x) — DN (u 2 )h{x)\ < ^ 11/3)1^ f |ui(y) - u 2 {y)\\h{x - y)\dy + ]- f \h(y)\\ui{x - y) - u 2 {x - y)\dy 

^ J Xq ^ J Xq 

PX l PX! 

+ \Kx )I Halloo / \My) ~My)\dy+\ui(x) -U 2 (x)\\\f \\ 00 | h(y)\dy 

J Xq j Xq 

An application of Young’s inequality for convolutions (see Theorem 2.24]) to the first two integrals gives, 

\DAf(ui)h(x) - DN{u 2 )h(x)\ < 11/3)1^ ||ui - u 2 || \\h\\ + \h(x)\ H/311^ ||ui - u 2 || + Mz) - u 2 (x )| ||/i|| . 

Consequently, taking the integral of both sides with respect to x yields 

\\DAf(ui)h(x) - DN(u 2 )h(x)\\ < (xi - x 0 ) H/311^ ||m - w 2 || \\h\\ + H/JH^ ||ui - u 2 || ||/i|| + ||«i - u 2 || H/3IU ||/i|| 
for all h € H. Then it follows that 


\\DJ\f{ui) - DJ\f{u 2 )\\ < (xi - ,T 0 + 2) H^ll^ Hm - u 2 || , 

which in turn implies that the nonlinear operator A f is continuously Frechet differentiable on H. □ 

For the zero solution, / = 0, we have DAf(f) = 0. Proposition [2] implies that in order to derive a sufficient condition, 
we also have to show that wo(£) < 0. To achieve that we follow the steps introduced in @, §6.1] and (To}. First, we show 
that Co semigroup (T(t)) t>0 generated by C is, in fact, a positive semigroup on the Banach lattice H with the usual sense 
of ordering “>” and the absolute value, | • |Q 

Proposition 4. The linear operator C on H = L 1 [to, Ti] generates a positive Co semigroup (T(t )) t>0 . Hence, the growth 
bound loo{C) is equal to the spectral bound s(£), i.e., s(£) = wo(£). 

Proof. For given A?C, the resolvent operator of C is given explicitly by 


Q\X) 

R( A, C)(f = (((x^gixo)-^ exp 


f 


A + w(s) 

ff(s) 


ds 


g{x, 


f X <p(y) exp (- f 

J Xq \ j y 


A + w(s) 
g(s) 


ds I dy 


(7) 


Using the similar reasoning to Lemma 2.4 of {H we conclude that for all A € C with Re( A) > + ||q|j , the resolvent 

operator R( A, C) exists. Thus, the resolvent set of C is not empty, i.e., p(C) ^ 0. Then, since g, w , and q are positive 
functions on [to, Ti], it is straightforward to see that R( A, C) is a positive operator whenever it exists. Consequently, 
the positivity of T(t) follows from the Characterization Theorem in [71, p.353]. Since H = L 1 [xo, x\\ is a Banach lattice 
and (T(t)) t > 0 is a positive semigroup, the last statement of the proposition follows from the main theorem presented in 

sm. □ 


1 We refer readers to [l3| for a detailed definition of a positive semigroup on a Banach lattice. 
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Recall that an operator A has a compact resolvent if p{A) ^ 0 and R( A, A) is compact for all A G p(A). Moreover, if 
one can prove that R( A, A) is compact for some A G p(A), then R( A, A) is compact for all A G p(A') . Operators which 
have a compact resolvent have the very nice property that their spectrum <r(A) coincide with their point spectrum a p (A), 
i.e., 

<j(A) = op(A) := {A G C | Acj) = A <f> for some </> ^ 0 G X} , (8) 

where X is some Banach space and A is an operator on X with a compact resolvent. 

Before proceeding further we first prove the following lemma, which we use in the proof of the next proposition. 

Lemma 5. For cf> G H and Ai G C with Re( Ai) > + HgH^ we define the linear operators 

rx 1 

Ki[tf>](x) = <j>{x 0 )g{x 0 ) = / q(y)<l>(y) dy 

J Xq 

and 

/v 2 [</>](£) = J <t>{y) exP f Xl d d ■ 

Then K\ and K 2 are compact on H. 


Proof. We first prove that K 2 is compact on H. Recall that an operator is compact if it maps bounded sets into relatively 
compact sets. In order to show that a set is relatively compact in L ?, we use the Kolmogorov-Riesz compactness theorem 
(see Ml for a statement of the theorem). 

Define a unit ball in H = L 1 [a;o, Xi] as B = {cj> G H | ||</>|| < 1}. Now, we prove that K 2 B is relatively compact by 
proving that each condition of the Kolmogorov-Riesz theorem applies to the set K 2 B. For a given function cf> G B, 


\K 2 <f>(x)\\ = [ [ 00) exp 

Jxo J X 0 


Ai + w(s) 


g{s) 


ds I dy 


dx < 


rx 1 rX 

/ / \4>{y)\ 

J Xn j Xc\ 


dy dx 


< (x\ - x 0 ) ll^ll < (xx - Xo) ■ 


This implies that the set K 2 B is bounded. So the first condition of the Kolmogorov-Riesz theorem is satisfied. Since we 
are working on a finite domain [xo, xi], the second condition of the Kolmogorov-Riesz theorem is satisfied by default. 

For convenience, let us define 


k(x, y) := exp 



g(s) J 


Observe that for Ai G C with f?e(Ai) > ||g|j + ||g|| , k(x, y) is uniformly continuous on [xo, x{\ x [xo, Xi]. So for a given 

£\ > 0 there exist (5i > 0 such that 


|fc(x + h, y) — k(x, y)\ < £\ for all \h\ < <5i and y G [xq, Xi] . 


(9) 


Furthermore, notice that \k(x, y)\ < 1 for x > y. Consequently, for a given e > 0 and <j> G B we have 


fX\ nX-\-h nX nX\ nX-\~n 

/ / k{x + h, y)<j>{y) dy- k(x, y)<t>{y) dy dx< / k(x + h, y)<t>{y) 

J Xn J Xn J xn J xc\ J x 


Xl nX-\-h 


dy 


Xq 

Xl 


' Xo 
0 X 1 


dx 


[k(x + h,y)- k(x, y)\ f>(y) dy 


dx 


< 


nx 1 px-\-h />xi px 

/ / \k(x + h, y)\ • \\<f>\\ dy dx + / £i\(j>(y)\dy dx 

J Xq J X J Xq J Xq 


< (Xl - X 0 )|/l| • || 0 || + £l(xi - X 0 ) • || 0 || 

< (^1 + £ 1)(x 1 - Xo) < £ 


for sufficiently small £\ and <5i in ©. This proves the third condition, and thus from the Kolmogorov-Riesz theorem the 
set K 2 B is relatively compact in H. This in turn implies that K 2 is a compact operator on H. 

For the operator Ki observe that the set I\\B is bounded, i.e., 


11-^100)11 




0(y)g(y) dy 


dx < (xi 


®o) Iklloo II011 < Oi -a:o) Iklloo ■ 


Since the function q(y), inside the integral, does not depend on x, the third condition of the Kolmogorov-Riesz theorem 
follows immediately. Hence, the operator K\ is also compact. □ 
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Proposition 6. The operator C defined in has a compact resolvent. 

Proof. From Proposition Q] we already know that the resolvent set of C is not empty. So we only have to prove that 
R( A, C) is compact for some A £ p(£). Particularly, we will prove that R( Ai, C) is compact for any Ai 6 C with 
.Re(Ai) > + Iklloo- The resolvent operator defined in ® can be written as the sum of the compositions of linear 

operators 

i?(Ai, C) = BiKi + B 2 K 2 , 

where K\ and K 2 are defined as in Lemma [5J 


Bi [4>\(x) 



g(s ) J 


4>{x ), 


and 

B 2 [<t>](x) = . 

g{x) 

From the assumptions (^) and (M,) it follows that the operators B\ and B 2 are bounded. In Lemma [5] we have proved 
that the operators K\ and K 2 are compact on H. Then the operator B\K\, composition of a bounded and a compact 
operator, is compact (see Proposition 5.43 of [HI). Similarly, the operator B 2 K 2 is compact. This in turn implies that 
the resolvent operator I?(Ai, £), a linear combination of compact operators, is also compact. □ 

The above proposition together with equation d8]) implies that the spectrum of the operator C consists of only eigen¬ 
values, i.e., cr(£) = a p (£). Thus, we can now characterize the spectrum of £ by its eigenvalues. 


Proposition 7. For A € C, 

where 

is a characteristic function of £. 


A £a p (£)=a(£) ^ £(A) = 0 , 


C(A) = [ 

J X 


g(x) 

xo g( x ) 


■ exp - 


f 

J X( 


A + w(s) 
g(s) 


ds 


dx — 1 


Proof. The eigenvalue equation for the operator C is given by 

Ccj) — A <f> = — ( g(x)cf>(x))' — w(x)<f>(x) — \<f> = 0. 


The solution of the above equation is given by the following eigenfunction 


(f{x) 


4>{xo)g(xo) 

g(x) 


exp 



A + w(s) 
g{s) 



( 10 ) 


We note that 


m = [ 

J X 

<-L 


Xq 

Xi 


\(j>(xo)g(x 0 )\ 

g{x) 

\<f(x 0 )g(x 0 )\ 


exp - 


exp 


ix 0 g{x) 

< {x\ - xo)\<f(xo)g(xo)\ 


-Re{ A) f 

J Xq 

r 1 

-Re A / — 

Jxn g{s) 


g(s) 


ds ) exp [ — 


/ 

J X( 


w(s) 

g(s) 


ds I dx 


ds I dx 


g(x) 


exp - 


n 1 

-Re( A) / — 

Jx 0 


ds 


and hence <f> £ H. From assumption {s^ g ) we have g £ C' 1 [:ro> x\]. This in turn implies that g’cf £ H. Analogously, we 
can also prove that (gf>)' = g'<f) + g<j>' £ H. 

Lastly, in order for (f £ V(C) we should have 


g(x 0 )(/)(x 0 ) = £[</>] = [ cj)(x 0 )g{x 0 )^-exp (- f A + w{s) \ 

Jx 0 g(x) V Jxo g{s) J 


which is equivalent to 




A + w(s) 

9(s) 



dx — 1 


□ 
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Since the eigenfunction tp defined in (flOl) is L 1 , the characteristic function £(A) is a finite-valued function for all A £ C. 
Moreover, when £(A) is restricted to R it is strictly decreasing. Furthermore, a simple limit calculation shows that 

lim £(A) = —1 and lim £(A) = oo . 

A—>oo A—)■—oo 


This in turn, from the Intermediate Value Theorem, implies that there exists a unique Ao £ R such that £(Ao) = 0. We 
can also guarantee that this eigenvalue Ao is negative real number provided that we have £(0) < 0. 

Remark 8. We also claim that the spectral bound of £ is equal to this Ao, i.e., s(£) = Ao- Suppose that there exists 
Ai £ er(£) = <T P {C) C C such that i?e(Ai) > A 0 . On the other hand, from Proposition [71 Ai should be a zero of 
characteristic function, i.e. £(Ai) = 0. However, 


1 = Re(£(Ai)) + 1 = f cos f 

Jx o 9i x ) YJx 


< 


< 


9(x] 

9(x) 

x 0 sO) 


XQ 

XI 


Im(Xi) 


ds 


L J x 0 
exp ( - 


exp - 


<?( 


x 


x 0 9(x] 

= 1 , 


exp - 


L 

L 


g(s) 

x Re{ Ai) + w(s) 
x 0 g{s) 

A 0 + w(s) 


f 


Ref Ai) + w(s) 


XO 

ds I dx 


g(s) 


ds dx 


9(s) 


ds dx 


which is a contradiction. 

We summarize the discussion of this section in the following criterion. 

Criterion. (Stability) The spectral bound of C is a unique real number Ao such that £(Ao) = 0 and hence ujo(jC) = s(£) = 
Ao • Moreover, the zero solution of the semilinear evolution equation © is locally exponentially stable if 


m = 



g(x) 

g(x] 


■ exp 



1 < 0 . 


3. Linearized instability for the zero solution 

In this section we derive sufficient conditions for instability of the zero solution. In particular, we use the following 
proposition from [13, p.206]. 

Proposition 9. Let (T(t)) t>0 be a Cq semigroup in the Banach space X with infinitesimal generator C. LetJf : X —> X be 
continuously Frechet differentiable on X. Let f £ 2A(£) be a stationary solution of ©. If there exists Ao £<7(£ + ZW(/)) 
such that Re(Xo) > 0 and 


max < oji(C + DAf(f)), 


sup 

\&a D (C+DM(f))\{\ 0 } 


Re{ A) > < Re(Xo) 


( 11 ) 


then f is an unstable equilibrium in the sense that there exists e > 0 and sequence {f n } in X such that f n f and 
\\ T (n)f n - f\\ > £ for n = 1 , 2 ,... . 

The discrete spectrum of an operator A on a Banach space X, denoted by <jd(A), is the subset of A £ a p {A) such that A 
is an isolated eigenvalue of finite multiplicity, i.e., the dimension of the set {if £ X : Aip = Xip} is finite and nonzero. Let 
(T(t)) t>0 be a Co semigroup on the Banach space A' with its generator A. Then the limit w\{A) = lim^oo f -1 log (a[T(t)]) 
is called a-growth bound of (T(i)) t > 0 . Here, a[T(t)] is a measure of non-compactness of T(t). The measure of non¬ 
compactness, introduced in a textbook [Ij|, associates numbers to operators (or sets), which tells how close is an operator 
(or a set) to a compact operator (or set). For instance, a[T(t)\ = 0 implies that the semigroup {T(t)) t>0 is eventually 
compact. In general, computing an explicit value of the a—growth bound uq is a complicated task. However, if we can 
prove that the linear operator C + DAf(f) generates an eventually compact Co semigroup, then from [2C| Remark 4.8] it 
follows that uq(£ + DJf(f)) = — oo. 


Proposition 10. The semigroup T(t ) generated by the operator £ is eventually compact. Specifically, it is compact for 
t > 2T(a;i), where 
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Proof. We first show that the semigroup ( T(t)) t>0 is differentiable for t > 2T(xi). Note that eventual differentiability 
implies eventual norm continuity (see Diagram 4.26 in [3, p. 119]). Consequently, the result follows from the compactness 
of the resolvent set of C (Proposition [Gj and [ 3 , Lemma 4.281. 

The following proof has been adopted from the proof of [9, Theorem 3.1]. The abstract Cauchy problem ut = Cu can 
be rewritten as a partial differential equation 


u t (t, x) + g(x)u x (t, x ) + (g'(x) + w{x)) u(t, x) = 0 . 

By Thoerem 2.2 in Q we know that the semigroup (T(t)) t>0 generated by C is given explicitly by 


T(t)<p(x) = 


Thus for t > T(x), we have 


ip (T 1 (T(x) - t)) t)] exp (- /r-i (r(x) _ t) da) t< T(x) 

^y/C (T (t - T(x)) <p(x)) exp (- /* o ds^j T(x) < t 

1 f Xl ( f x w(s ) \ 

i(t, x) = -r-rr / q(x)u(t — T(x), x)dx x exp — / ; ds ] . 

9 { x ) J xo \ J xo g(s) ) 


( 12 ) 


(13) 


On the other hand, since g € C 1 [xo, x\] and g(x ) > 0 we have T(x) € C\x o,xi]. Therefore, u(t, x) is continuous both in x 
and t for t > T(xi), where T(xi) is the maximum of T(x) on [xo, Xi]. Then, from (fl2l) . u(t, x) is continuously differentiable 
for t > 2T(xi). □ 

In Proposition [6] we have proved that the spectrum of C consists of only eigenvalues, which can be expressed as the 
zeros of the characteristic equation £(A). Consequently, similar to Remark [8] we can show that there exist a unique real 
number Ao > 0 such that 


sup Re{ A) < sup Re( A) < Re{ A 0 ) 

Aea D (£+DA^(/))\{Ao} Ae,r p (£+DA^(/))\{A 0 } 


if and only if 


m = f 

J X. 


X1 q{x) 


exp 


7(S) 


/xo 9i x ) \ Jx 0 g(s) 

We can now summarize the results of this section in the following criterion. 


ds) — 1 > 0 . 


Criterion. (Instability) The a—growth bound of the operator C can be found explicitly, i.e., ui\(C) = — oo. Furthermore, 


/ xq \ j xo 1 ?(^) / 

then the zero solution of the semilinear evolution equation ([3]) is unstable in the sense described in Proposition [3 


4. Concluding remarks 

Expressing stability results in terms of characteristic functions has been popular since A. Lotka published his pioneering 
article on population modeling 17]. The characteristic function, £(A), that we derived in this paper is easy to compute. For 
example, consider a hypothetical scenario for which the model parameters are given by xo = 1, Xi = 1000, q(x) = ln(x), 
g{x) = Yg.x(1001 — x), w(x) = Y(jgo(£ — l) 1 ' 17 . In this case the zero solution of the evolution equation ([1]) is locally 
asymptotically stable. On the other hand, decreasing the growth rate twofold gives an unstable zero solution. This might 
not be very intuitive at first glance. However, when we decrease the growth rate twofold, aggregates will grow slower, but 
new cells keep entering single cell population at the same rate (fecundity rate q{x)). Thus, to keep the zero solution stable 
we should also decrease the fecundity rate (at least) twofold. Furthermore, as one might expect, doubling the fecundity 
rate q[x) also makes the zero solution unstable. 

As a future research, we plan to extend the results of this paper to nontrivial solutions (stationary, self-similar, etc.) 
of this aggregation-growth model as well as apply this linearization method to a size-specific aggregation-fragmentation 
model considered in 16]. 
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